Chandra Observations of the Central Region of Abell 3112 



Motokazu Takizawa 1 ' 2 , Craig L. Sarazin 1 , Elizabeth L. Blanton 1,3 , and Greg B. Taylor 4 



O 



ABSTRACT 



We present the results of a Chandra observation of the central region of Abell 



3112. This cluster has a powerful radio source in the center and was believed to 
have a strong cooling flow. The X-ray image shows that the intracluster medium 
(ICM) is distributed smoothly on large scales, but has significant deviations from 
a simple concentric elliptical isophotal model near the center. Regions of excess 
emission appear to surround two lobe-like radio-emitting regions. This structure 
probably indicates that hot X-ray gas and radio lobes are interacting. From an 
1 analysis of the X-ray spectra in annuli, we found clear evidence for a temperature 

decrease and abundance increase toward the center. The X-ray spectrum of the 
central region is consistent with a single-temperature thermal plasma model. 
^ ■ The contribution of X-ray emission from a multiphase cooling flow component 

with gas cooling to very low temperatures locally is limited to less than 10% 
of the total emission. However, the whole cluster spectrum indicates that the 
ICM is cooling significantly as a whole, but in only a limited temperature range 
(> 2 keV). Inside the cooling radius, the conduction timescales based on the 
Spitzer conductivity are shorter than the cooling timescales. We detect an X-ray 
point source in the cluster center which is coincident with the optical nucleus of 
v~j ■ the central cD galaxy and the core of the associated radio source. The X-ray 

spectrum of the central point source can be fit by a 1.3 keV thermal plasma 
and a power-law component whose photon index is 1.9. The thermal component 
is probably plasma associated with the cD galaxy. We attribute the power-law 
component to the central AGN. 
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1. Introduction 

The central regions of clusters of galaxies are very interesting and enigmatic places. The 
intracluster medium (ICM) in the central region is often so dense that its radiative cooling 
time is significantly shorter than the Hubble time. Therefore, the ICM will cool down unless 
it is heated significantly. Unbalanced radiative cooling would cause a "cooling flow" in 
the cluster center (see Fabian 1994 for a review). X-ray imaging analyses by ROSAT and 
EXOSAT indicated that mass deposition rates, M, were more than 100M o yr~ 1 for many 
clusters (Edge, Stewart, & Fabian 1992; Allen & Fabian 1997; Peres et al. 1998). However, 
the fate of the resulting cold gas, or whether ICM in "cooling flow" clusters really cools down 
to low temperatures, is still unclear. The cD galaxies in cooling flows have an excess blue star 
component, which suggests that the cooled ICM forms new stars (McNamara & O'Connell 
1989). However, the implied star formation rates are much lower than the cooling flow mass 
deposition ones (McNamara & O'Connell 1989). Furthermore, recent XMM-Newton high 
resolution spectroscopic observations have imposed a strong constraint on the amount of 
gas cooling down to low temperatures (Peterson et al. 2001; Kaastra et al. 2001; Tamura 
et al. 2001), which is consistent with the fact that the mass deposition rates determined 
spectroscopically by ASCA tended to be lower than those determined by EXOSAT and 
ROSAT (see Makishima et al. 2001 for a review). 

The cD galaxies in the cooling flow clusters usually host strong radio sources, with 
a central radio core and/or radio lobes. Observations utilizing Chandra's excellent spatial 
resolution have revealed interactions between these radio lobes and ICM, which may affect 
the dynamical and thermal history of ICM in the cluster center. In Hydra A (McNamara 
et al. 2000; David et al. 2001), Perseus (Fabian et al. 2000), and Abell 2052 (Blanton et al. 
2001), radio lobes coincide with X-ray low surface brightness regions, which suggests that the 
radio plasma expands to displace the ICM. On the other hand, X-ray cavities without radio 
emission are also found at larger distances in Abell 2597 (McNamara et al. 2001). They are 
believed to be old radio lobes which cannot emit (at least at high radio frequencies) because 
the high energy relativistic electrons have already suffered radiative losses. These radio 
quiet cavities also indicate that cD galaxies provided the ICM with high energy electrons 
and magnetic fields intermittently. 
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Abell 3112 is a cooling flow cluster at a redshift of z — 0.0746. There is a powerful 
radio source, PKS 0316-444, in the cluster center. Previous X-ray imaging analyses with 
EXOSAT (Edge, Stewart, & Fabian 1992) and ROSAT (Allen & Fabian 1997; Peres et al. 
1998) indicated that Abell 3112 had a strong cooling flow, with a mass deposition rate of 
M ~ 400 M & yr _1 . The temperature profile at a larger scale than the cooling radius was 
studied by Markevitch et al. (1998) and by Irwin, Bregman, & Evrard (1999) with ASCA 
data. Markevitch et al. (1998) found that the temperature decreased outward from 6 keV to 
3 keV with ASCA data, although they explicitly assumed a cooling flow component in the 
central region. On the other hand, Irwin et al. (1999) argued that the data were consistent 
with an isothermal distribution. An increase in the iron abundance toward the center was 
reported by Finoguenov, David, & Ponman (2000). However, there was no detailed study 
of the temperature and abundance structure within the cooling radius (< 250 kpc) because 
of the limited spatial resolution of ASCA. Furthermore, any interaction region between the 
radio source and the ICM would be too small to have been resolved with ROSAT or ASCA. 
In this paper, we present Chandra observations of central region of Abell 3112. We assume 
H = 50 km s" 1 Mpc -1 , fi = 1-0, and fi A = 0.0. At a redshift z = 0.0746, 1" corresponds 
to 1.94 kpc. The errors correspond to the 90% confidence level throughout the paper. 



2. Observations and Data Reduction 

2.1. X-ray Observations 

Abell 3112 was observed twice with Chandra on 2001 May 24 for 7,257 seconds and on 
2001 September 15 for 17,496 seconds. Each observation was pointed so that the cluster 
center would fall near the aimpoint of the ACIS-S3 detector. However, the roll angles were 
different, so the outer parts of the field of view differ for the two observations. The two 
observations were merged, using the positions of bright point sources to register the two 
images. Here, we analyze data from the S3 chip only. The data from the SI chip was 
used to search for the background flares. A few short periods with background flares were 
found and removed using the LC_CLEAN software provided by Maxim Markevitch 5 , leaving 
a total exposure of 21,723 seconds. Only events with ASCA grades 6 of 0, 2, 3, 4, and 6 were 
included. Standard bad pixels and columns were removed. Exposure maps and background 
data were generated for each observation separately, and then merged appropriately. Blank 
sky background data were taken from the compilation by Maxim Markevitch 5 . 



5 See http://asc.harvard.edu/cal/Links/Acis/acis/Cal.prods/bkgrnd/current/index.html. 
6 See http://asc.harvard.edu/ciao/. 
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The data are possibly affected by the low energy QE degradation of ACIS. 7 We tried 
to correct this using the CORRARF 7 program. However, the correction did not appear to 
work well for this data; specifically CORRARF seemed to over-correct the data. As a result, 
we obtained absorption column densities which were much lower than the Galactic value 
in our spectral fits. We expect that the absorbing column should be at least as large as 
the Galactic value. Therefore, we will show the results without the correction for the low 
energy QE degradation. Fortunately, this correction significantly affects only the very low 
energy (< 1 keV) part of the Chandra band. As a result, neither temperature nor abundance 
determinations were seriously affected by this correction. On the other hand, mass deposition 
rates of cooling flow components increased by a factor of about two when we fitted the spectra 
after the correction. 



2.2. Radio observations 

The observations were made with Very Large Array 8 at a center frequency of 1320 MHz 
on 1996 October 18. The array was in the 'A' configuration, which provided an angular 
resolution of 6.9 x 1.5 arcsec. A total of 1.2 hours were obtained on source using 32 channels 
across a 6.25 MHz band. Both right and left circular polarizations were observed. Phase 
and bandpass calibration was obtained by short (1 min) observations of the strong (4.9 Jy) 
calibrator J0440-4333 taken every ~20 minutes. Atomic hydrogen absorption was searched 
for towards the compact core over a velocity range from 22,100-23,440 km s _1 at a resolution 
of 44 km s -1 , but not detected down to a 3cr limit on the maximum optical depth of 0.012. 



3. X-ray Images 

The X-ray image was adaptively smoothed using the ciao 9 csmooth program with 
a minimum signal-to-noise of 3 per smoothing beam. An identically-smoothed background 
was subtracted, and the image was divided by an identically-smoothed exposure map. The 
exposure map was corrected for vignetting, using the typical energy of the observed photons 
to take account of the energy dependence of the vignetting. The resulting image is shown 

7 See http://cxc.harvard.edu/cal/Links/Acis/acis/Cal_prods/qeDeg/index.htnil . 

8 The National Radio Astronomy Observatory is operated by Associated Universities, Inc., under cooper- 
ative agreement with the National Science Foundation. 

9 See http://asc.harvard.edu/ciao/. 
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in Figure 1. As in other typical cooling flow clusters, X-ray emission is distributed symmet- 
rically and is fairly concentrated in the center. The image is smooth and quite symmetric 
on large scales. The cluster emission is elliptical, with a major axis along the direction from 
NNE to SSW. On large scales, there appears to be no substructure except for some point 
sources. 

We used a wavelet detection algorithm (WAVEDETECT in CIAO) to detect individual 
point sources. Sources were visually confirmed on the X-ray image, and a few low-level 
detections were removed. Using this method, we found sixteen point sources. One of these 
corresponds to the radio source PKS 0316-444, which is located at the center of the central cD 
galaxy ESO 248- G 006. Another corresponds to a galaxy in the cluster (LCRS B031619.9- 
442441). 

The regularity of the cluster X-ray emission on large scales doesn't continue into the 
very central region. Figure 2 shows the X-ray contour map of the ~ 1' x 1' region around the 
central point source. In the very center, an elongated structure from the central point source 
towards the SE is clearly seen. We can see another filament-like component from ~ 10" 
south of the central point source towards the SW. To clarify structures which deviate from 
symmetric cluster emission, we fitted the entire cluster image with a concentric elliptical 
isophotal model, and then subtracted the best-fit model from the image to get the residual 
component. In the fitting procedure, we fixed the center of the isophotes to the center 
determined from the larger scale isophotes about which the image was symmetric, rather 
than the central point source. The ellipticity and position angle for each isophote were 
allowed to vary. It is most likely that relaxed cluster potential structure is not spheroidal 
but triaxial (Jing & Suto 2001). In this case, the ICM hydrostatic density distribution also 
will be triaxial (Lee & Suto 2003). Therefore, the ellipticities and the position angles of 
X-ray isophotes can vary with radius unless our line-of-sight coincides with one of the three 
symmetry axes. 

In Figure 3, the greyscale shows the residual produced by subtracting our best-fit el- 
liptical isophotal model from the adaptively smoothed image of the central region of the 
cluster. The sharp elliptical edge near the outside is an artifact because the elliptical model 
extended only to this isophote. Dark areas are positive residuals (excess emission), while 
light areas are negative residuals. The dark and light regions at the center are due to the 
central point source, which is not exactly at the center of the X-ray isophotes at large radii. 
There are regions of excess emission to the south of the central source. The contours show 
the 1.32 GHz VLA radio image. There is a very bright radio core, which is unresolved in 
the radio image; the elliptical contours for the core show the beam of the radio observation. 
The very bright radio core is coincident with the central X-ray source. There are two diffuse 
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radio regions to the SE and SW of the radio core, which may be connected to the core. 
The NS extent of these regions is presumably exaggerated due to the elongated observing 
beam. Roughly speaking, the excess X-ray regions appear to surround two lobe-like radio 
components. However, part of the SW radio lobe is coincident with a region of X-ray ex- 
cess emission. Unfortunately, the low resolution of the radio data, bright radio core, and 
elongated radio beam make it difficult to decide definitely whether the excess is in the radio 
lobes or partly surrounding it. 

4. Spectral Analysis 

In order to examine the temperature and abundance structure quantitatively, we de- 
termined the X-ray spectra in annuli. We fitted the data between 0.5 and 10.0 keV with 
a photo-absorbed single temperature MEKAL model (Kaastra 1992; Liedahl, Osterheld, & 
Goldstein 1995) after masking the point sources. We tried to fit the data fixing the absorbing 
column density to the Galactic value (N H = 2.51 x 10 20 cm~ 2 ; Dickey & Lockman 1990) 
or allowing it to vary freely. The fitting results depend somewhat on whether we fixed the 
absorption to the Galactic value or not. Figure 4 shows the radial temperature and abun- 
dance profiles. The solid crosses are the values obtained when we allowed the absorption to 
vary, while the dashed crosses are the values when we fixed the absorption to the Galactic 
value. The dashed crosses are shifted slightly to the left in order to be more easily seen. The 
results with fixed or varying absorption are similar. At r > 70", both the temperature and 
abundance are nearly constant at kT ~ 6 keV and Z ~ 0.4 Z e , respectively. The tempera- 
ture decreases from ~ 6 keV at 70" to 3.5 keV in the central 20". The abundance increases 
from ~ 0.4 Z Q at 70" to 1.3 Z Q in the central 20". 

The temperature values obtained with variable absorption, along with their associated 
errors, are fitted to a power-law profile given by T(r) = T (r/l") p . For r < 70", T = 
1.838 ± 0.254 keV and p = 0.270 ± 0.044. For r > 70", T = 4.226 ± 3.174 keV and 
p = 0.064 ± 0.165. We will use these fitted temperature profiles later when we analyze 
conduction timescale and integrated gravitational mass. 

To examine the non-radial temperature structure, we made a two-dimensional temper- 
ature map of the central 240" x 240" region. We divided the region into 64 (8 x 8) square 
regions, each of which is 30" x 30". Then, we fitted the data of each region with the photo- 
absorbed MEKAL model where absorption is treated as a free parameter. Figure 5 shows a 
two-dimensional temperature map of the central 240" x 240" region overlaid with X-ray sur- 
face brightness contours. Black and white represent lower (~ 3 keV) and higher (~ 6 keV) 
temperatures, respectively. There is no significant azimuthal structure in the temperature 
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map while temperature decrease toward the center is again clearly seen. We also made a 
two-dimensional abundance map in the same way. However, abundances in the individual 
fits were too poorly constrained to provide a useful map. 

We tried to fit the annular spectra with a model for a multiphase cooling flow to constrain 
the contribution from a cooling flow component and a mass deposition rate in each annulus, 
although the data within each annulus are almost consistent with a single temperature 
MEKAL model. As a cooling flow spectral model, we used the MKCFLOW model based on 
the MEKAL model. We added a MEKAL model representing the emission from the ICM 
outside the cooling flow. Both MKCFLOW and MEKAL were assumed to be subject to the 
same absorption column, which we fixed to the Galactic value or let vary freely. We fixed 
the metallicity and initial gas temperature of the MKCFLOW model to the same values of 
metallicity and temperature as the MEKAL model, respectively. The low temperature in the 
MKCFLOW model is fixed to the lowest value (0.001 keV) in XSPEC. The results are shown 
in Table 1 and 2. Compared with the fitting results without cooling flows, the results are not 
improved dramatically. Figure 6 shows a radial profile of the mass deposition rate derived 
from the fitting. The circles and solid error bars are the values obtained from fits with freely 
varying absorption. The squares and dashed error bars are the values from the fits with fixed 
Galactic absorption. When absorption is set as a free parameter, forth bin and the outermost 
three bins are consistent with a model without a cooling flow component. Only the inner 
three bins and the fifth bin have the cooling flow component, although their errors are fairly 
large and comparable to themselves. The resultant total mass deposition rate is AA^H'l M Q 
yr _1 , which is significantly smaller than the rates derived from the previous image analyses 
of EXOSAT and ROSAT data, which are ~ 400 M© yr" 1 (Edge, Stewart, & Fabian 1992; 
Allen & Fabian 1997; Peres, et al. 1998). The cooling flow component contributes less 
than ~ 10% of the total X-ray emission. When we adopt the Galactic absorption value, the 
mass deposition rate becomes even lower and consistent with a model without a cooling flow 
component in all bins. 

Even though multiphase gas is not found within small regions of the cluster, we do see an 
overall temperature gradient, with gas in the center that is cooler by a factor of ~ 2 compared 
with gas in the outer regions of the cluster. Therefore, if we fit a spectrum extracted from 
the whole cluster, a significant mass deposition rate and a low temperature cutoff might be 
expected. We fit the data within 157" with the MEKAL plus MKCFLOW model where the 
low temperature in MKCFLOW is allowed to vary. The results are shown in Table 3. The 
upper and lower rows show the results with the absorption column allowed to vary and fixed 
to the Galactic value, respectively. In contrast to the spectral analyses in annuli, the mass 
deposition rate is comparable with the former values based on imaging analysis. However, 
the low temperature in the MKCFLOW model is not very low (/cTl ow ~ 2 keV), which 
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means that the observations are consistent with significant cooling, but only over a limited 
temperature range. This is consistent with what has been observed with many cooling flow 
clusters using Chandra and XMM- Newton (e.g., Peterson et al. 2003). Please note that the 
goodness of the fit is marginal (the reduced chi-squared is about 1.5). This might indicate 
that the emission measure distribution with the temperature is not that of the standard 
cooling flow model even in the temperature range between Tmgh and Tl ow - 



To quantify the radial structure of the ICM, we made a radial profile of the X-ray surface 
brightness in the 0.3-10.0 keV band (Figure 7a). The bins for X-ray surface brightness were 
chosen to be the same as those used below to determine the radial variation in the X-ray 
spectrum (§ 4). The surface brightness values were deprojected to determine the X-ray 
emissivity and gas density (Figure 7b), assuming the emissivity is constant in spherical 
shells. When we convert the emissivity into the gas density, we use a MEKAL code and 
assume that the temperature is constant in spherical shells and equal to that derived from 
the spectral fitting of projected data. Some fluctuations are seen in the density and pressure 
profiles near the outer boundary (Figure 7b and 7c), which are artifacts of the deprojection 
since we assume zero emissivity outside the outer boundary. However, many simulations 
of our deprojection method show that this affects only a few outermost points, and has no 
significant effect anywhere near the center. The observed gas density and pressure profiles 
are quite smooth outside of the region of the radio source. Based on the smooth radial 
profiles and lack of substructure in the X-ray image outside of the radio source region, we 
conclude that this cluster is dynamically relaxed and that the ICM is generally in hydrostatic 
equilibrium. 



The gravitational mass distribution was determined from the equation of hydrostatic 
equilibrium, 



We use the fitted power-law temperature profiles when M tot (r) is calculated. To reduce 
the noise in the density gradient term, we determine the gradient at by differencing the 
densities at rj_2 and r^. Uncertainties in the density gradient at are also calculated from 
the density uncertainties at rj_ 2 and r i+2 . The density gradients and their uncertainties are 



5. X-ray Profiles and Deprojection 



6. 



Mass Distributions 




(1) 
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assumed to be equal in the innermost three points and the outermost three points. The gas 
mass profile is also determined from the ICM density profile assuming spherical symmetry. 
The gas mass is then given by 

M gas (<r) = / p gas (r)4nr 2 dr (2) 
Jo 

The integrated gravitational mass and gas mass are shown in Figure 8 (a). 

A radial profile of the gas mass fraction is shown in Figure 8 (b). The gas mass fraction 
increases from 3% at 10" to 10% at 100". This trend is similar to those in other well-relaxed 
clusters (David, Jones, & Forman 1995; Ettori, & Fabian 1999; Allen, Schmidt, & Fabian 
2002). Our gas mass fraction is lower than that of Mohr, Mathiesen, & Evrard (1999) at 
much larger radii (> 1 Mpc), but the gas fractions of clusters generally increase with radius, 
and this appears to apply to Abell 3112. Using the same beta-model fits as in Mohr et al., 
we find that the gas fraction at 100" is predicted to be ~12%, which agrees well with our 
values. 



7. Comparison of Cooling and Thermal Conduction 

Because the temperature is higher than ~ 3 keV over the whole cluster, the main emis- 
sion mechanism is thermal bremsstrahlung. Therefore, the isobaric cooling time is (Sarazin 
1986) 

Radiative cooling would make a radial temperature gradient in the central region of the 
ICM through its density dependence, as is observed. (Fig. 4). On the other hand, thermal 
conduction would have the effect of reducing the temperature gradient. The conduction 
timescale is generally expressed as (see also Sarazin 1986) 

tcond = ^ (4) 
K 

where n e is electron number density, It = T/VT is the scale length of the temperature 
gradient, and ks is the Boltzmann constant. If we consider only the Coulomb scattering 
process, the thermal conductivity for hydrogen plasma is (Spitzer 1962) 

K = 46 x 1 ° u (wk) " 2 (ir) erg " cm " K " ■ (5) 

where In A, the Coulomb logarithm, is 



In A = 37.8 + In 



(l0 8 K:)(lO- 3 cm- 3 ) .' ^ 
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Figure 9 shows radial profiles of t coo i and t con d- When we calculated t con d, we used the 
fitted power-law temperature profile to reduce the noise in the temperature gradient. Also, 
tcond is evaluated only for r < 70" because an isothermal distribution is consistent with the 
temperature data in the outer regions. 

The cooling time of the innermost bin is 2 x 10 9 yr. If we define a cooling radius as 
where t coo \ is 1.5 x 10 10 yr, it becomes ~ 100", or ~ 200 kpc. Both results are consistent 
with the former ROSAT (Allen & Fabian 1997; Peres et al. 1998) and EXOSAT results 
(Edge, Stewart, & Fabian 1992). At the innermost three bins, t coo i and t con d are comparable 
with each other. Outside of this region, the conduction timescale is clearly shorter than the 
cooling timescale. If conduction proceeded at the Spitzer rate, it would erase the temperature 
gradient. Since the observations show a significant gradient in this region (Fig. 4a), thermal 
conduction must be significantly suppressed by, e.g., tangled magnetic fields and/or plasma 
instabilities. The cooling rate of gas with solar abundance is only ~ 10% higher than that 
of hydrogen-helium mixture gas at 3.5 keV (see Figure 9-9 of Binney & Tremaine 1987). 
Therefore, actual cooling timescales in the central regions are probably ~ 10% shorter at the 
most, which does not change the situation here dramatically. Although our estimation here 
is rough, it is still useful for an estimation of an order of magnitude. Calculations taking 
account of time evolution will be helpful to make more precise models. 



8. The Central Point Source 

We detected an X-ray point source in the cluster center. The position is coincident with 
the optical core of the cD galaxy and the radio core. The position of the X-ray point source 
is (epoch J2000) a = 03 /l 17 m 57.65 s 5 = -44°14'17.06" with an uncertainty of 0.04 arcsec in 
each coordinate. The spectrum of the point source is shown in Figure 10. We fitted the data 
with a MEKAL model plus a power-law model. In addition to the absorption column density 
common for both components, we added intrinsic absorption for the power-law component 
only. As a result, the applied model is 

Model AG N = WABS x [MEKAL + ZABS x POWERLAW]. (7) 

The abundance of the MEKAL component was fixed to the value derived from the region 
surrounding the point source. The absorbing column A% was also fixed either to the value 
derived from the surrounding region or the Galactic value. The fitting results are shown 
in Table 4. The spectral model consists of a 1.26 keV thermal plasma component and a 
power-law component with a photon index of ~ 1.9. The power-law component has an extra 
absorbing column of ~ 2 x 10 21 cm -2 . This is consistent with the central source being a 
strongly absorbed AGN. 



9. Discussion 



We found that the central radio source associated with the central cD galaxy in Abell 3112 
is interacting with the surrounding ICM (Fig. 2 and 3). The radio image shows a central core 
and two diffuse lobes to the SE and SW of the core. The central X-ray image has an asym- 
metric structure, and the excess X-ray emission over an elliptical isophotal model roughly 
appears to surround two radio lobes. This configuration is naturally explained if the radio 
lobes have swept up the ICM which was where the radio lobes are now. Similar scenarios 
have been proposed in other cooling flow clusters with radio sources (e.g., McNamara et al. 
2000; Fabian et al. 2000; Blanton et al. 2001 ). In addition, there is excess X-ray emission 
coincident with part of the SW radio lobe. It is possible that the excess X-ray emission is 
due to the cool gas that was originally near the very center; this cool gas might have been 
entrained by hot buoyant bubbles. This scenario was originally proposed to explain the X- 
ray morphology of the central region of M87 (Churazov et al. 2001). Indeed, the filamentary 
structures in the X-ray image of the central region of Abell 3112 (Figure 2) are similar to 
those in M87 (Feigelson et al. 1987; Bohringer et al. 1995). This model may also apply 
to Abell 133 (Fujita et al. 2002). In order to determine which model applies to Abell 3112, 
temperature and abundance measurements of the excess X-ray components would be useful 
because the excess components are expected to have higher abundances and lower entropy 
in the entrainment scenario. Unfortunately, the present data have too few counts in these 
regions to allow a detailed spectral analysis. A deeper radio image would also be useful, as 
would lower frequency observations to search for steep spectrum emission. Note that the 
effect of the central radio source is limited to a region very close (r ~ 10") to the center. This 
region is much smaller than the cooling flow region itself (r ~ 100"). Outside the interacting 
region, there is no evidence of the interaction such as X-ray cavities, instead the X-ray image 
is quite symmetric and smooth. Therefore, the central radio source does not directly affect 
the structure of the cooling flow outside of the very central region, at least at present. Note 
that we only have high frequency radio data. Low frequency radio observations might give 
another picture. In the case of M87 in the Virgo cluster, for instance, 327 MHz radio data 
show an morphology which fills out a large portion of the central region of the cluster (Owen, 
Eilek, & Kassim 2000), although similar structures cannot be seen in higher frequency data. 

Radiative cooling produces a temperature decrease towards the center of the cluster. 
On the other hand, thermal conduction would reduce this temperature gradient. If we adopt 
the standard Spitzer conductivity, conduction would be expected to be very effective and 
to have eliminated the central temperature gradient (Figure 9). However, the observations 
show that the temperature does decrease significantly at the center of the cluster. This 
implies that thermal conduction is significantly suppressed below the Spitzer value. 
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Although we did not find a large amount of locally cooling multiphase gas, the existence 
of a large-scale temperature gradient indicates that the ICM is cooling in a global scale. The 
mass deposition rate of Abell 3112 determined from our spectroscopy of the total spectrum 
is comparable to that determined previously based on R OS A T and EXOSAT image analyses 
(Edge, Stewart, & Fabian 1992; Allen & Fabian 1997; Peres et al. 1998). However, the 
temperature range of the cooled gas is limited to > 2 keV. In addition, emission measure 
distribution about temperature is probably different than that expected from the standard 
cooling flow model. This suggests that some other heating mechanism may affect the cooling 
gas. One probable solution is heat conduction which is reduced from the Spitzer value by 
some physical mechanism but still is energetically important (e.g., Hattori & Umetsu 2000; 
Malyshkin & Kulsrud 2001). The relativistic electrons from the central radio source might 
also supply the required heating (e.g., Fabian et al. 2000; McNamara et al. 2000; David et 
al. 2001; Blanton et al. 2001). However, it is unlikely that the radio-emitting electrons are 
heating all of the cooling gas in Abell 3112 at present, because the radio lobes are so small that 
they do not directly affect the cooling flow on large scales. In Abell 3112, there is no evidence 
for ghost bubbles which might indicate a history of past radio source interactions (McNamara 
et al. 2001). Another possible form of radio source heating would be high energy protons 
from the central AGN. Since the minimum energy density in relativistic electrons is too small 
to support radio bubbles against the surrounding hot gas (e.g., Blanton et al. 2001), most 
of the energy in the radio source may be in the form of relativistic protons. The accelerated 
protons could diffuse and heat the ICM through Coulomb interactions (Rephaeli & Silk 1995; 
Inoue & Sasaki 2001). Although they do not emit any observable radio radiation, protons 
would produce several hundred MeV 7-rays through the decay of neutral pions produced in 
collisions with thermal protons. Future 7-ray observations will test this hypothesis. 
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Fig. 1. — Adaptively smoothed Chandra image (0.3-10.0 keV) with a minimum signal-to- 
noise of 3 per smoothing beam. The image has been background subtracted and corrected 
for exposure and vignetting. 
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Fig. 2. — X-ray contour map (0.3-10.0 keV) of the central ~ 50" x 70" region of the cluster 
Abell 3112. Contours are logarithmically spaced. In the very center, an elongated structure 
is seen extending from the central point source towards the SE. Another filament extends 
from ~ 10" south of the central point source towards the SW. 
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Fig. 3. — The greyscale shows the residuals of the adaptively-smoothed X-ray image (0.3-10 
keV) from our best-fit elliptical isophotal model. Dark areas are positive residuals (excess 
emission), while light areas are negative residuals. The dark and light regions at the center 
are due to the central point source, which is not exactly at the center of the X-ray isophotes 
at large radii. The contours are from the 1.32 GHz VLA radio image. The very bright radio 
core is coincident with the central X-ray source. The NS extension of the radio core source 
is due to the elongated radio observing beam. 
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Fig. 4. — Radial profiles of temperature (a) and abundance (b), respectively. The solid 
crosses are the values obtained when we allowed the absorption to vary, while the dashed 
crosses are the values when we fixed the absorption to the Galactic value. 
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Fig. 5. — A two-dimensional temperature map of the central 240" x 240" region overlaid with 
X-ray surface brightness contours. Black and white represent lower (~ 3keV) and higher 
(~ 6keV) temperature, respectively. 
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Fig. 6. — Radial profile of the mass deposition rate. The circles and solid error bars are the 
values obtained when we allowed the absorption to vary, while the squares dashed error bars 
are the values when we fixed the absorption to the Galactic value. 
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Fig. 7. — Radial profiles of surface brightness (a), electron density (b), and pressure (c), 
respectively. 
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Fig. 8. — (a) Integrated total gravitational mass (diamonds) and gas mass (triangles), (b) 
Radial profile of gas mass fraction. The gas mass fraction increases from 3% at 10" up to 
10% at 100". 
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Fig. 9. — Radial profiles of isobaric cooling timescales (circles) and thermal conduction 
timescales (squares). At the innermost three bins, cooling and conduction timescales are 
comparable with each other. Outside this region, conduction timescales are clearly shorter 
than cooling timescales. 
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Fig. 10. — Spectrum of the central point source. The spectrum is fit by a 1.26 keV thermal 
plasma plus a power-law component whose photon index is 1.9. The power-law component 
has an intrinsic extra absorption column of ~ 2 x 10 21 cm" 2 . This is consistent with the 
central source being a strongly absorbed AGN. 
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Table 1: Cooling Flow Spectra in Each Annulus with Variable Absorption 
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Table 2: Cooling Flow Spectra in Each Annulus with Fixed Absorption 
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Table 3: Total Spectrum (r < 157") 







Z 


iV H 


M 


X 2 /dof 
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(z & ) 


(10 20 cm- 2 ) 
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Note. - The upper and lower rows show the results with absorption column density allowed 
to vary and fixed to the Galactic value, respectively. 
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Table 4: Spectrum of the Central Point Source 



kT 


z 


N n 
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x 2 /dof 
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(Ze) 


(10 20 cm- 2 ) 


(10 20 cm" 2 ) 




1.2618:3 


(0.77) 


(6.40) 


1 SQ+0-30 1 8 Q+18.1 


12.9/21 




(1.27) 


(2.51) 


2.00^ 25.3+ 1 f 


12.9/21 



Note. - The abundance and Nh of the MEKAL component were derived from the fit of the 
region surrounding the point source. Upper and lower rows show the results when the iV H is 
allowed to vary or fixed to the Galactic value, respectively. 



